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ABSTRACT 


In  Section  I  an  attempt  is  made  to  derive  the  degeneracy  factor 
of  Eyring's  Significant  Structure  Theory  for  a  particular  model  of  a 
solid  with  vacant  lattice  sites.  A  form  is  obtained  which  differs 
from  that  obtained  by  Eyring  and  which  gives  results  less  satisfactory 
than  Eyring's  factor  when  compared  with  Hugoniot  data. 

In  Section  II  some  equilibrium  Hugoniot  states  are  calculated  for 
a  mixture  of  two  materials  and  the  formulae  are  applied  to  combinations 
of  quartz-like  and  polyethylene-like  materials.  Gruneisen  parameters 
are  found  to  differ  somewhat  from  those  calculated  with  a  simple  mass- 
weighted  mixing  formula  and,  in  general,  the  properties  of  the  more 
compressible  material  tend  to  dominate.  Some  elastic  calculations  are 
made  for  various  composite  configurations.  It  is  noted  that  conditions 
of  uniaxial  strain  are  violated,  and  that  the  relations  between  mean 
pressure  and  volume  are  sensitive  to  the  boundary  condition  on  the  face 
of  the  composite.  Some  comments  are  made  about  deviations  from  equili¬ 
brium. 


Conventional  methods  for  reducing  ultrasonic  velocity  data  taken 
at  high  pressures  require  a  priori  knowledge  of  the  compressibility  as 
a  function  of  pressure  if  elastic  constants  are  to  be  calculated.  In 
Section  III  an  iterative  procedure  is  outlined  which  makes  it  possible 
to  determine  the  elastic  constants  provided  only  that  ambient  pressure 
values  of  thermal  expansion  coefficient,  its  temperature  derivative,  and 
specific  heat  at  constant  pressure  are  known.  The  method  is  applied  to 
travel  time  measurements  in  NaCl  and  KC1  made  by  Bartels  and  Schuele. 

The  new  procedure  changes  values  of  the  elastic  constants  by  about  1%. 


Section  IV  extends  the  above  iterative  procedure  to  include  data 
and  constants  for  crystals  of  arbitrary  symmetry.  It  is  applied  to 
cal  cite  and  rutite  for  illustration. 
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PART  I 


EQUATION  OF  STATE  FROM  SIGNIFICANT  STRUCTURE  THEORY 

C.  T.  Tung 


Two  computations  have  been  made  in  attempts  to  improve  the  Signifi¬ 
cant  Structure  Theory  described  in  Technical  Summary  Report  No.  3.  The 
theory  requires  at  every  P  and  T  the  computation  of  the  "solid"  part  of 
the  volume,  V$,  before  the  liquid  volume,  V,  can  be  determined.  In  Sum¬ 
mary  Report  No.  3,  V$  is  computed  from  an  extrapolation  of  the  experimental 
solidus  for  argon.  It  appears  logically  superior  to  use  the  solid  part 
of  the  partition  function  to  compute  Vs,  even  though  this  is  out  of  the 
equilibrium  range  of  the  solid  phase.  The  program  was  modified  to  accomp¬ 
lish  this,  but  the  results  were  unsatisfactory,  apparently  because  the 
solid  partition  function  is  unsatisfactory  for  outside  the  equilibrium 
region. 

In  the  Significant  Structure  Theory  of  liquids,  a  partiton  function 
of  a  liquid  is  considered  as  three  parts:  solid  partition  function,  degen¬ 
eracy  factor,  and  gas  partition  function.  Two  proportionality  constants  are 
introduced  into  the  degeneracy  factor  of  Eyring's  model  and  they  are  deter¬ 
mined  by  fitting  the  thermodynamic  properties  at  normal  conditions;  it  is 
doubtful  that  these  two  constants  are  really  constant  over  a  wide  range  of 
pressures  and  temperatures.  To  avoid  this  confusion  and  to  search  for  a 
more  realistic  mathematical  description  of  the  degeneracy  factor,  we  give 
an  elementary  derivation  of  this  factor  for  simple  liquids  (such  as  argon, 
xenon,  neon,  krypton,  etc.)  based  on  some  simple  geometrical  configurations. 

Let  us  define  the  following: 

<j> ( r )  —  intermolecular  potential  of  separation  distance  r. 

V„  —  molar  volume  of  solid, 
s 

V  —  molar  volume  of  liquid. 

X  =  Vs/V  —  probability  for  a  site  to  be  occupied  by  a  molecule. 
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Y  =  (V-V  )/V  --  probability  for  a  site  being  empty  (hole). 
M  --  coordination  number  (=12). 
a  --  nearest  neighbor  distance. 


As  far  as  a  given  atom  A  is  concerned,  the  positional  degeneracy  of 
atom  A  is  essentially  determined  by  its  nearest  neighboring  populations  of 
atoms  and  holes.  There  are  thirteen  types  of  geometrical  configurations  as 
follows  (actual  picture  being  three-dimensional): 

(•  denoting  atom  and  o  denoting  hole) 


type  j  = 

0 

1 

j 

12 

Ba 

Y~p\ 

W 

1  A  J 

V  A  J 

vlv 

configuration 

no  hole 
adjacent 
to  A 

1  hole 
adjacent 
to  A 

j  holes 
adjacent 
to  A 

12  holes 
adjacent 
to  A 

pj 

1 

1 

probability 
of  jth  type 

c«xmy° 

cmxm-iyi 

i 

cV_JYJ 

J 

c"x°ym 

ej 

total  energy 
of  atom  A  in 
jth  type 

Eo. 

E-|  =  EQ-4>(a) 

Ej  =  E0-J4>(a) 

Em  =  E0-M*(a) 

Eq  is  the  total  energy  of  atom  A  when  the  first  shell  is  completely  filled. 
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Let  us  denote  the  averaged  total  energy  of  atom  A  over  twelve  con¬ 
figurations  as  E.  Then,  we  have 

E  =  1  PiE1  =  1  c|2X12-j[E  -  j 4> ( a ) ]Y J 

j=0  J  J  j=0  J 

=  E  ]f  C1.2X12_jYj  -  4>(a)  l  jcl2X12'jYJ 
j=0  J  j=l  J 


=  EQ(X  +  Y)12  -  <f>(a)Y[9(X  +  Y)12/aY] 


=  Eq  -  12Y<f>(a)  =  Eq  -  12*(a)(V  -  Vs)/V  (1.1) 

Now,  let  us  observe 

exp(-E/kT)  =  [exp(-E0/kT)][exp(+12<{)(a)(V  -  Vg )/VkT) ] . 

Thus,  we  get 

positional  degeneracy  factor  =  expC+lZ^a) (V  -  V$)/VkT].  (1-2) 

If  |+12<f>(a) (V  -  V$)/VkT|  «  1 ,  we  have 

positional  degeneracy  factor  =  1  +  1 2<f>(a) ( V  -  Vs)/VkT.  (1.3) 

Here,  nearest  neighbor  distance,  a,  can  be  expressed  in  terms  of  Vg  by 

a  =  (v^Vs/N)^3,  N  being  Avogadro's  number. 

The  validity  of  Eq.  (1.2)  or  Eq.  (1.3)  is  restricted  to  the  condition  that 
<j>(a)  be  positive,  since  the  positional  degeneracy  factor  must  always  be 
greater  than  or  equal  to  one. 
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Replacing  Eyring's  degeneracy  factor,  g  -  1  +  n^exp[-aEsVs/(V-V$)NkT] , 

by  Eq.  (1.2),  we  calculate  the  Hugoniot.  The  results  are  less  satisfactory 
than  those  obtained  with  the  Eyring  model,  indicating  that  his  factor  sum¬ 
marizes  effects  other  than  the  spatial  degeneracy  considered  here. 


PART  II 


COMPRESSION  OF  COMPOSITE  MATERIALS 
G.  E.  Duvall  and  S.  M.  Taylor,  Jr. 


2.1  INTRODUCTION 

We  define  a  composite  material  as  one  in  which  the  physical  properties 
vary  from  point  to  point,  usually,  but  not  necessarily,  periodically.  A 
lattice  of  mass  points  with  interacting  forces  is  an  example  of  such  a  sys¬ 
tem,  and  the  study  of  wave  propagation  in  lattices  is  very  old.  Sir  Isaac 
Newton  studied  the  propagation  of  waves  in  a  one-dimensional  lattice  as 
substitute  for  a  continuous  string.  With  mathematical  tools  available  to 
him  he  found  the  lumped  constant  system  to  be  simpler  than  the  continuum. 

Sir  William  Hamilton  also  studied  the  one-dimensional  lattice  and  found 
the  equivalent  of  the  Schroedinger  solution,  though  it  was  not  so  concisely 
expressed  since  Bessel  functions  had  not  yet  been  invented.  In  the  latter 
part  of  the  nineteenth  century  Lord  Baden-Powell  and  later  Lord  Kelvin  used 
the  one-dimensional  lattice  as  a  model  for  the  study  of  optical  dispersion. 

In  modern  physical  science  and  technology,  problems  of  x-ray  diffrac¬ 
tion,  electronic,  mechanical  and  thermal  properties  of  solids,  E-M  antennae 
design,  and  even  astrophysics  all  hinge  upon  the  effects  of  wave  propagation 
in  periodic  or  aperiodic  structures  or  systems. 

Problems  in  the  present  context  vary  somewhat  from  those  conventionally 
associated  with  lattices  because  it  is  not  always  possible  to  ignore  details 
of  propagation  within  the  components  of  the  composite  material  and  because 
problems  of  finite  amplitude  waves  are  often  involved.  This  means  that  the 
dispersion  phenomena  which  are  associated  with  any  lattice  become  more  com¬ 
plicated  and  that  the  superposition  of  multiple  disturbances  is  no  longer 
permitted.  In  spectral  terms  this  means  that  propagation  velocity  depends 
on  both  frequency  and  amplitude  and  that  waves  of  a  single  arbitrary  frequency 
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cannot  be  propagated.  In  such  circumstances  it  is  sometimes  helpful  to 
consider  and  attempt  to  understand  a  very  specific  and  well-defined  simple 
problem,  which  may  then  serve  as  a  base  for  the  understanding  of  more  com¬ 
plicated  problems. 

For  this  purpose  we  consider  a  half-space  of  the  composite  material 
and  suppose  that  the  free  surface  is  subjected  to  a  step-change  in  pressure 
or  velocity  normal  to  the  surface.  This  state  of  pressure  or  velocity  at 
the  free  surface  is  maintained  indefinitely  and  we  may  inquire  about  the 
development  and  structure  of  the  wavefront,  deviations  from  conditions  of 
uniaxial  strain,  and  the  equilibrium  states  in  a  region  long  after  the  wave- 
front  has  passed. 

Simple  as  this  problem  is,  its  complete  solution  is  still  very  ambi¬ 
tious  and  we  confine  ourselves  for  the  present  to  questions  about  the 
equilibrium  state,  with  some  remarks  about  the  developing  wavefront. 


2.2  EQUILIBRIUM  COMPRESSION  OF  TWO  FLUIDS 

A.  B.  Wood  many  years  ago  considered  a  composite  system  consisting  of 

2  1 

small  gas  bubbles  dispersed  in  water.  *  He  argued  that  at  low  frequencies 
there  should  be  local  equilibrium  of  temperature  and  pressure  amongst  gas 
bubbles  and  water,  that  each  should  assume  a  volume  appropriate  to  that 
pressure  and  temperature,  and  that  the  derivative  of  the  resulting  pres¬ 
sure-volume  relation  would  accordingly  give  the  sound  velocity.  To  a  first 
approximation  the  water  could  be  assumed  rigid,  then  the  sound  velocity  is 


c  *  (P/p0R)1/2(1  +  R) 


(2.1) 


where  P  is  ambient  pressure,  pQ  =  density  of  water  and  R  is  the  volume  ratio 
of  gas  to  water  at  P  and  temperature  T. 


As  would  be  expected,  sound  velocities  predicted  by  Eq.  (2.1)  can  be 
very  low,  as  little  as  a  few  meters/second.  Equation  (2.1)  has  a  broad 
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minimum  at  R  =  1  with  c  =  10  meters/sec  when  P  =  1  bar,  with  c  rising 

steeply  for  both  large  and  small  R  (Fig.  2.1).  H.  B.  KarpTus,  in  work 

reported  about  ten  years  ago,  refined  Wood's  model  and  tested  it  experi- 

2  2  2  1 

mentally  for  considerable  range  of  parameters.  '  ’  *  His  results  are 
summarized  in  Fig.  2.1,  taken  from  reference  2.2.  The  frequencies  at  which 
measurements  were  made  were  less  than  10%  of  the  mean  resonant  frequency  of 
the  bubbles,  and  he  concluded  that  there  was  no  significant  frequency  dis¬ 
persion. 

The  principles  on  which  Wood's  model  are  based  are  well  founded  and 
can  be  extended  to  other  materials.  When  a  compression  front  passes  over  a 
mixture  of  two  materials  there  will  be  an  initial  period  in  which  stress 
equilibrium  is  being  established.  The  duration  of  this  period  will  be  the 
order  of  several  times  the  larger  of  the  two  numbers  a./c^,  i  =  1,  2,  where 
sl  and  c  are  characteristic  dimension  and  sound  velocity  of  the  two  compo¬ 
nents  respectively.  As  a  result  of  compression,  each  component  will 
experience  a  temperature  rise  determined  by  the  stress  field  and  by  its 
own  material  properties.  Roughly  speaking,  when  stress  equilibrium  has  been 
established,  temperatures  T-j  and  will  exist  in  the  two  components.  Neg¬ 
lecting  convection  and  radiation,  which  will  accelerate  the  process,  these 

two  temperatures  will  be  equalized  in  a  time  the  order  of  several  times  the 

2 

larger  of  the  two  numbers  C.p^Ji.  /k^  ,  i  =1,  2,  where. C.»  are 

specific  heat,  density,  dimension  and  thermal  conductivity  of  the  ith  com¬ 
ponent  respectively.  The  time  required  for  thermal  equilibration  is 
normally  much  greater  than  that  required  for  stress  equilibration.  The 
ratio  of  the  latter  time  to  the  former  is  the  order  of  i/h,  where  a  is  the 
scattering  mean  free  path  for  phonons.  A  may  range  from  a  single  lattice 
cell  diameter  in  a  disordered  structure  to  grain  size  in  a  highly  ordered 
structure. 

Sufficiently  far  behind  the  front  of  the  shock,  then,  we  expect  stress 
and  temperature  equilibrium  to  exist.  If  the  materials  be  liquid,  this 
reduces  to  equilibrium  in  P  and  T.  In  either  liquid  or  solid  case  the 
thermodynamics  of  compression  are  appropriately  described  by  the  Gibbs  func¬ 
tions  of  the  two  components  and  of  the  mixture.  If  interfacial  energies  are 
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VELOCITY  OF  SOUND, M/S 


VELOCITY  OF  SOUND 
IN  WATER  CONTAINING  AIR  BUBBLES 

Figure  2.1 
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negligible,  we  can  disregard,  for  the  present,  the  sizes  and  shapes  of  the 
components  and  consider  only  their  mass  fractions.  Let  X  be  the  mass  frac¬ 
tion  of  component  2  and  1  -  x  that  of  component  1.  Then  the  Gibbs  function 
for  unit  mass  of  the  mixture  is 

G  =  (1  -  X)  G1  +  X  G2  (2.2) 
where  subscripts  "1"  and  "2"  refer  to  the  two  components.  Other  thermody¬ 


namic  functions  are  obtained  as  derivatives  of  Eq.  (2.2): 

specific  volume:  V  =  (9G/9P)T  =  Gp  (2.3) 

=  (1  -  X)V1  +  XV2  (2.4) 

specific  entropy:  S  =  -(9G/9T)p  =  -Gy  (2.5) 

=  (1  -  x)S7  +  XS2  (2.6) 

isothermal  compressibility: 

3  =  “(1/V) ( 9 V/ 9  P ) T  =  -Gpp/V  (2.7) 

ev  =  (l  *  x)31v1  +  xe2v2  (2.8) 

thermal  expansion  coefficient: 

a  =  (l/V)(9V/9T)p  =  GpT/V  (2.9) 

aV  =  (1  -  X)a1V1  +  Xa2V2  (2.10) 

specific  heat  at  constant  pressure: 

Cp  =  T(9S/9T)p  =  -TGtt  (2.11) 

Cp  =  (1  -  x)Cpl  +  xCp2  (2.12) 
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G  and  its  derivatives  are  assumed  to  be  continuous,  so  the  definitions  in 
Eqs .  (2.7),  (2.9)  and  (2.11)  lead  to  the  following  identities,  which  can 
be  of  use  in  evaluating  experimental  data: 

GppT  =  -(9f5V/aT)p  =  (3aV/3P)T  (2.13) 

GTTP  =  -0/T)(3Cp/3P)T  »  (3(aV)/3T)p:  (2.14) 

In  the  application  of  Eqs.  (2.2)  -  (2.12),  we  assume  the  thermodynamic  pro¬ 
perties  of  the  two  components  to  be  known.  These  equations  provide  procedures 
for  computing  properties  of  the  mixture. 

If  the  compression  front  envisioned  in  this  problem  ever  achieves  what 
can  be  described  as  a  steady  profile,  then  the  jump  conditions  can  be  applied 
to  relate  the  equilibrium  state  behind  the  front  to  that  ahead.  The  thermo¬ 
dynamics  of  compression  are  then  summarized  in  the  Rankine-Hugoniot  Equation, 

E  =  E0  +  (1/2) (P  +  PQ)(V0  -  V),  (2.15) 

and  in  the  first  and  second  laws  of  thermodynamics.  Differentiating  Eq. 
(2.15),  we  have 

dE  =  (1/2) (VQ  -  V)dP  -  (1/2) (P  +  PQ)dV  =  TdS  -  PdV 


or 

TdS  =  (1/2) (VQ  -  V)dP  +  (1/2) (P  -  PQ)dV  (2.16) 

Expressing  both  V  and  S  in  terms  of  P  and  T,  as  in  Eqs.  (2.3)  and  (2.4)  we 
obtain  for  the  equation  of  the  Hugoniot  in  the  T,  P  plane: 

dT  =  dP[(VQ  -  V)  -  3V(P  -  P0)  +  TaV]/[2Cp  -  aV(P  -  PQ)]  (2.17) 

Integration  of  Eq.  (2.17)  makes  it  possible  to  evaluate  any  thermodynamic 
parameter  in  the  equilibrium  compressed  state,  provided  G^(P,T),  G2(P,T) 
and  A  are  known. 
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The  Gruneisen  parameter  r  and  sound  velocity  c  are  readily  obtained 
from  the  above  expressions.  From  the  definition  of  sound  velocity: 

c2  =  -V2(3P/9V)s  =  -V2/[(aV/aP)s],  (2.18) 

(9V/aP)s  =  (aV/aT)p  +  V(aT/aP)s 

=  -eV  +  Ta2V2/Cp  (2.19) 

Substituting  (2.19)  into  (2.18)  yields 

c2  =  V2/(bV  -  Ta2V2/Cp)  (2.20) 

The  Gruneisen  parameter  is  calculated  in  similar  fashion: 

r  =  (V/cv)(aP/aT)v  =  Va/cve  (2.21) 

Since  Cy  =  Cp  -  a2VT/g, 

Eq.  (2.21)  becomes 

r  =  aV2/(eVCp  -  a2V2T).  (2.22) 

In  computing  c  and  r  for  the  mixture,  V,  aV,  gV  and  Cp  are  to  be  determined 
from  Eqs.  (2.4),  (2.8),  (2.10)  and  (2.12),  then  substituted  into  Eqs.  (2.20) 
and  (2.22). 

In  illustrating  the  above  relations  and  procedures,  we  would  like  to 
say  that  we  had  created  Gibbs  functions  for  two  substances  which  were  con¬ 
sistent  with  all  available  experimental  data  and  had  all  of  the  proper 
limiting  behavior.  Unfortunately  we  can't.  What  we  have  done  is  to  create 
rather  conventional  descriptions  of  a  quartz-like  material  and  a  polyethylene¬ 
like  material  and  use  these  to  illustrate  the  effects  of  varying  X  on  the 
above  parameters.  The  following  assumptions  were  made  for  each  material: 
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(2.23) 

(2.24) 

(2.25) 

(2.26) 


Polyethylene 

Quartz 

i 

1 

2 

b.j ,  g/cc 

.3887 

.7832 

Cvi ,  Mbcc/g° 

1.436  x  10"5 

.56  x  10 

Vio»  cc/9 

1.035 

.378 

Yio 

10.2 

5.58 

19.31 

2.70 

The  results  of  computations  are  shown  in  Figs.  2.2  -  2.6.  Most  are 
not  particularly  remarkable.  Figure  2.2  shows  that  temperature  along  the 
Hugoniot  increases  very  rapidly  when  a  little  polyethylene  is  added  to 
quartz,  and  that  the  rate  of  change  of  temperature  decreases  as  the  amount 
of  polyethylene  increases.  This  is  the  kind  of  effect  to  be  expected  when 
a  compressible  material  is  added  to  a  relatively  incompressible  one.  The 
nearly  linear  change  in  V/VQ  for  fixed  P  shown  in  Fig.  2.3  shows  that  vol¬ 
ume  is  but  little  affected  by  temperature  changes  on  the  Hugoniot,  again  as 
expected.  In  Fig.  2.4  the  Gruneisen  parameter,  computed  from  Eq.  (2.22), 
is  compared  for  =  0.50  with  a  simple  mass-weighted  average: 

r  =  (1  -  x)r1  +  xr2  (2.27) 
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Figure  2.2 


Temperatures  and  Pressures  on  the 
Hugoniots  of  Mixtures  of  Quartz 
and  Polyethylene. 

X  =  mass  fraction  of  quartz 
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/x  / 

X<7  5 
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P,  Megabars 


Figure  2.4 


Figure  2.5  shows  perhaps  the  most  interesting  result.  Although  r/V  is 
assumed  constant  for  each  component,  the  resultant  for  the  mixture  is  quite 
strongly  volume  dependent;  the  mass  weighted  average  values  differ  appre¬ 
ciably  from  the  correct  values.  Sound  velocities,  shown  in  Fig.  2.6,  change 
almost  linearly  with  X. 

These  results  on  variations  of  r  and  r/V  with  volume  and  with  x  are 
worthy  of  particular  note.  Attempts  have  been  made  elsewhere  to  characterize 
multicomponent  materials  in  terms  of  effective  r's  in  such  a  way  that  the 
effective  r  turns  out  to  be  approximately  equal  to  or  less  than  that  of 
the  lower  constituent.  Experiments  performed  to  check  these  calculations 
have  produced  pressures  lower  than  calculated  values  by  as  much  as  56  per¬ 
cent,  indicating  that  the  effective  Gruneisen  parameter  was  even  lower  than 
estimated.  This  result  is  contrary  to  those  shown  in  Fig.  2.4;  however, 
experiments  were  performed  with  slightly  porous  samples,  and  this  fact  calls 
the  results  into  question. 

The  influence  of  varying  r  on  the  stress  field  has  been  the  subject 
of  extensive  debate.  It  seems  fairly  certain  that  in  some  circumstances  it 
can  have  a  substantial  effect.  For  that  reason  it  is  most  appropriately  cal¬ 
culated  by  the  procedure  given  here  when  equilibrium  conditions  are  expected. 
If  they  are  not,  it  is  of  no  significance  anyway. 

2.3  LOCAL  DEVIATIONS  FROM  UNIAXIAL  STRAIN 

The  feature  of  plane  stress  wave  propagation  in  homogeneous  materials 
which  most  clearly  distinguishes  it  from  other  geometries  is  the  exact  main¬ 
tenance  of  uniaxial  strain;  no  material  particle  of  a  continuum  undergoes 
any  motion  in  a  direction  parallel  to  the  wavefront.  This  is  one  of  the 
first  conditions  to  fail  in  a  composite  material,  except  in  a  special  case 
discussed  below. 

To  illustrate  the  point  we  consider  several  special  cases: 
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Figure  2.6 


2.3.1  Two  Fluids;  Slabs  in  Propagation  Direction 

The  fluid  structure  is  assumed  to  be  symmetric  so  that  on  com¬ 
pression  from  xQ  to  x  the  cell  width,  aQ,  remains  constant.  The  total 
specific  volumes  before  and  after  compression  are  given  by  Eq.  (2.4).  Apply¬ 
ing  this  to  the  system  of  Fig.  2.7  yields  the  relations 

Vo  =  V2(P-T)X 

(a0  "  wQ)x0  =  V7 (P,T) (1  -  x) 

wx  =  V2(P  +  dP,  T  +  dT)x 

(a  -  w)x  =  Vy(P  +  dP,  T  +  dT) (1  -  x) 

If  we  choose  xQ  so  the  mass  of  the  cell  is  one  gram,  let  x  -  x  =  dx,  w  -  wQ 
=  dw,  V2(P  +  dP,  T  +  dT)  -  V^P.T)  =  dV2,  etc.,  the  lateral  strain  resulting 
from  compression  dx  is 

(2.28) 
(2.29) 


e  2  =  dw/wQ  =  dV2/V2  -  dV/V 
=  -dw/(aQ  -  wQ)  =  dV1/V]  -  dV/V 


These  satisfy  the  cell  condition  f  0  e  dy  =  0. 

0  V 


2.3.2  Two  Fluids;  Slabs  Parallel  to  Wave  Front 

In  this  case,  shown  in  Fig.  2.8,  lateral  strain,  ey,  is  iden- 

|  '  tically  zero,  dx  =  dV,  dw  =  dV?. 

L 


2.3.3  Two  Fluids;  Rods  Normal  to  Wave  Front 

Symmetry  produces  a  rectangular  cell,  of  dimensions  yQ  by  zQ, 
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Compression  of  Two  Fluid  Slabs 
Normal  to  Wave  Front 

/  . 

Figure  2.7 
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Compression  of  Two  Fluid  Slabs 
Parallel  to  Wave  Front 


Figure  2.8 

i 
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on  the  boundary  of  which  the  particle  motion  is  normal  to  the  plane  of  the 
paper  (Fig.  2.9).  Again  choosing  xQ  so  as  to  produce  a  cell  of  unit  mass, 
we  have 


x  v  z  =  V 
0  0  0 


Tra  X  =  A  x  =  xV0 

oo  oo  2 


(yozo  •  Ao>xo  =  (1  •  x)vl 


After  compression  x  =  xQ  +  dx,  a  =  aQ  +  da,  A  =  AQ  +  dA. 


dx  =  dv/y0z0 


A„dx  +  x  dA  =  xdV0 
o  o  2 


(y0z0  '  A0)dx  "  xodA  =  ^  '  x)dV] 


The  relative  change  in  area  of  the  rod  is 


dA/A0  =  (dV2/V2  -  dV/V) 


(2.30) 


If  it  remains  circular  in  cross-section,  the  radial  strain  is 


da/a0  =  dA/2A0 


(2.31) 


2.3.4  Two  Fluids;  Rods  Parallel  to  Wave  Front 


The  mass  per  unit  length  of  the  cell  shown  in  Fig.  2.10  is 


M  =  x0y0/V. 


That  of  the  cylinder  is 


MX  = 


VV2 


"a0/V2 
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Compression  of  Rods 
Normal  to  Wave  Front 

Figure  2.9 
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Compression  of  Rods 
Parallel  to  Wave  Front 


Figure  2.10 


When  pressure  P  is  applied  as  shown,  the  area  of  the  rod  changes  by  an 
amount  dA  =  MAdV2;  if  it  retains  its  circular  cross-section  the  radial 
strain  at  its  periphery  is 

da/aQ  =  dV2/2V2  (2.32) 


The  same  result  is  obtained  if  the  rods  are  staggered  so  that 
the  undistorted  cell  is  hexagonal. 


2.3.5  Elastic  Slabs  Normal  to  Wave  Front^’^ 

In  this  case,  shown  in  Fig.  2.11,  it  is  not  enough  to  apply 
a  uniform  pressure  in  the  direction  of  compression.  Instead  we  suppose  the 
cell  to  be  compressed  by  a  rigid  piston  driven  by  a  force  F  per  unit  length. 
Then  the  mean  pressure  acting  over  the  cell  is  p"  =  F/aQ.  Under  the  assumed 
conditions  the  x-displacement  is  the  same  for  both  components,  so 

Exl  =  ex2  ■  dx/xo 

In  the  notation  of  Fig.  2.11  we  define  specific  volumes  V-| ,  V2,  V: 

V  •  a0x0;  (1  -  x)V,  =  (a0  -  w0)x0;  xVj  =  v.Qx0 

Differentiating  these,  we  have 

aQdx  =  dV;  wQdx  +  xQdw  =  XdV2*,  (aQ  -  wQ)dx  -  xQdw  =  (1  -  X )dV^ 

For  small  dx: 


ey2  5  dw/W0  -  dVV2  -  Ex;  Ex  "  dx/X0 

eyl  '  'dw/(ao  -  wo>  =  dVl/Vl  ‘  Ex 
ezl  £z2  =  0 
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Compression  of  Elastic  Slabs 
Normal  to  Wave  Front 


Figure  2.11 


Assuming  Hooke's  law  with  Lame  constants  L-j ,  y-| ,  L2» 

°xl  =  Lld'fl/Vl  +  2"1ex 
cj  I  =  (L-j  +  2vi-|)dV-|/V-|  -  2^1 


°2l  * 

°x2  =  ^2d^2^2  +  2u2ex 

ay2  =  ^*"2  +  "  2y2ex 

oz2  =  L2dV2/V2 


In  addition  we  have  the  relations 


dV  =  XdV2  +  (1  -  x)dV] 

°yi  =  V 

0xl^  "  A^V1  +  ax2xV2  =  "Fxo 

Define 

|/1  =  (1  -  X )  ;  l/2  =  XV2;  a  =  L  +  2y 

Q  =  +  a2/l/2 

Then  the  above  equations  yield  the  following: 


dVj/Vj  =  ex[l  +  (L2  -  L-|)/V-|Q] 
dVV2  =  ex[1  '  (L2  _  Ll)/l/2Q] 

VI  =  'x<L2-L1^V 


(2.33) 

(2.34) 

(2.35) 
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ey2  =  ■£x^L2  '  L1^2^ 

(2.36) 

°xl  =  *  L1<L?  " 

(2.37) 

ayl  =  ex^Ll  +  “ 

(2.38) 

a2l  =  li£x[i  +  a2  -  l1>/v1q: 

(2.39) 

ax2  =  ex^a2  +  L2^L1  "  L2^/V2Q^ 

(2.40) 

°y2  =  ex[L2  +  a2^Ll  " 

(2.41) 

az2  =  W1  +  <L1  -  L2>/%« 

(2.42) 

To  a  first  approximation,  waves  will  be  transmitted  down  the 
cities  given  by 

slab  with  velo- 

C-|  -  V-]9oxi/8ex  -  [ct-j  +  L-j  (L2  ~  L-|)/(/-|Q]Vi 

(2.43) 

C2  =  V23ax2/,9ex  =  ^-a2  +  L2^L1  "  L2^I/2^V2 

(2.44) 

For  comparison,  the  velocity  of  waves  on  a  thin  plate  in  air 

is  given  by 

Cp  =  (a  -  U)/p 

(2.45) 

where 

Poisson's  ratio  =  5/(1  +  5) 

£  plays  the  same  role  in  extension  of  a  plate  with  ez  =  0  as  Poisson's  ratio 
in  extension  of  a  bar.  Comparing  Eqs.  (2.43)  and  (2.44)  with  (2.45)  shows 
that  effective  values  of  c  can  be  defined  for  each  component: 

51  eff  =  (L1  ‘  L2^/l/lQ 

(2.46) 

?2  eff  =  ^L2  " 

(2.47) 
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If,  for  example,  >  L-| ,  ^  <  0  and  component  ”1”  contracts  laterally 

when  F  >  0.  As  -»•  0,  eff  ?2  =  L2^a2’  ancl  comP°nent  "2"  becomes  an 
unsupported  plate.  When  l_2  =  L-| ,  uniaxial  strain  exists. 

If  components  2  and  1  have  large  and  small  moduli,  respec¬ 
tively,  compression  in  the  x-direction  produces  a  resolved  shear  stress, 
ax2  "  ay2’  1n '2  which  will  ultimately  lead  to  structural  failure  of  2.  A 
non-zero  modulus  in  component  1  delays  failure  of  2  when  loading  is  increased. 
For  example,  let  failure  occur  when 

<J>  =  (cx  -  cry )2  +  (ox  -  az)2  +  (0^  -  az)2  =  2k2  (2.48) 

Then  for  an  unsupported  plate  with  given  e 

A 

*p  =  +  S2  +  o  +  C)2].  (2.49) 

When  5  =  Keff>  $  =  4>ef f ;  when  £  =  0»  uniaxial  strain  obtains  and  g  =  0, 

4>  =  «(>0 .  Then  when  L2  > 


*0  -  *2  eff  -  V 

Failure  then  occurs  for  a  value  of  a  less  than  the  failure  value  in  uni- 

A 

axial  strain  and  greater  than  that  for  an  unsupported  plate.  Under  the 
same  conditions,  5-j  <  0  and  <f>^  <  <j>0,  so  the  value  of  ax  at  which 

failure  of  1  occurs  is  greater  than  for  uniaxial  strain. 

The  net  result  of  all  this  is  that  waves  in  the  softer  compo¬ 
nent  will  travel  faster  than  in  uniaxial  strain,  those  in  the  stiffer 
component  will  travel  more  slowly  than  in  uniaxial  strain,  and  that  there 
will,  in  general,  be  two  cusps  in  the  (p  ,V)  or  (a.e  )  curve  for  the  com- 

A  A  A 

posite  due  to  elastic  failure  of  the  components. 

The  above  discussion  has  been  cast  in  terms  of  linear  elastic 
materials  for  convenience  and  because  the  results  are  explicit.  The  results 
are  readily  generalized  if  e  and  a  are  replaced  by  de  and  da,  respectively. 
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2.4  DEVIATIONS  FROM  EQUILIBRIUM 

As  stated  in  Section  2.2,  thermal  equilibration  takes  longer  than 
stress  equilibration.  The  time  required  depends  upon  dimensions  and  mate¬ 
rials;  and  in  many  real  and  experimental  situations,  thermal  equilibrium 
will  not  be  reached.  Fortunately  wave  propagation  is  sometimes  insensitive 
to  thermal  effects,  unless  they  produce  phase  transitions  or  chemical  reac¬ 
tions.  In  any  exceptional  cases  which  must  be  treated,  the  effects  can  be 
bounded  by  assuming  isothermal  compression  on  the  one  hand  and  adiabatic 
on  the  other. 

For  inert,  fixed  inhomogeneities,  stress  equilibration  time  ^  <x£/c, 

1  <  a  <  10,  will  usually  be  a  reliable  guide.  If  the  inhomogeneities  are 
not  inert  or  are  not  fixed,  a  new  set  of  rules  applies.  The  most  notable 
example  is  in  the  effect  of  dislocations  on  precursor  decay.  The  mechanics 
of  any  such  case  must  be  examined  in  detail  in  order  to  make  any  estimate 
of  equilibration  time.  Another  practical  example  of  such  a  case  is  that  of 
a  granular  material  in  which  small  grains  may  slide  over  one  another  under 
the  influence  of  shear  forces  and  friction. 

There  is  another  manner  in  which  the  usual  concepts  of  equilibration, 
discussed  in  Sections  2.2  and  2.3,  may  be  violated:  In  Fig.  2.11  an  ele¬ 
ment  of  a  cell  is  assumed  to  be  compressed  uniformly  by  a  rigid  piston  so 
that  ex  is  the  same  for  each  component.  If,  in  fact,  a  real  x-travelling 
shock  were  being  produced,  one  component  would  be  compressed  more  than  the 
other,  resulting  in  relative  motion  at  the  interfaces  and  in  different  stress 
conditions  from  those  obtained  in  Section  2.3.5,  unless  the  shock  were  too 
weak  to  break  the  interface  bond. 

Consider,  for  example,  alternate  slabs  of  quartz  and  polyethylene  in 
the  geometry  of  Fig.  2.12.  A  uniform  pressure,  p-j,  is  applied  to  the  top 
of  the  column  and  held  there.  Wave  propagation  velocity  in  quartz  is  greater 
than  in  polyethylene,  so  waves  in  the  center  slab  outrun  those  on  the  sides. 

A  true  steady  state  is  never  reached.  However  a  kind  of  equilibrium  may  be 
achieved.  The  compression  wave  causes  the  quartz  slab  to  expand,  radiating 
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Figure  2.12 
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energy  into  the  polyethylene  ahead  of  the  polyethylene  wave  front.  This 
tends  to  inhibit  the  rate  of  propagation  of  energy  in  quartz.  If  we  sup¬ 
pose  that  this  effect  is  great  enough,  we  may  assume  the  jump  conditions 
to  be  approximately  satisfied  across  a  wavefront  travelling  down  the  cell 
with  velocity  D  in  both  components.  Then  the  jump  conditions  for  mass  and 
momentum  become 

°oi“oiD  "  Vi*0  -  ui>  (2'50) 

Vi  =  ”oi*o1Dui  (2'51) 

1  =  ’-2i  *ol  -  ao  '  V  *02  =  wo 

*1  =  ao  '  *2  =  " 

Define  =  (aQ  -  w)x-j ,  =  wx2,  where  x-|  and  x2  are  the  respective 

lengths  in  "1"  and  "2"  of  elements  which  originally  had  lengths  xQ.  Define 

xrt  so  that  the  mass  of  the  undisturbed  cell  of  volume  x  a  is  unity.  Then 
o  oo 

the  jump  conditions  can  be  summarized  in  the  relation 


v^/v2  =  0  -  X)x1(x0  -  X])/  x2(x0  -  x2)  (2.52) 


If  we  examine  the  stress  and  strain  in  the  compressed  region,  we  get  a  set 
of  relations  identical  to  those  set  forth  in  Section  2.3.5,  except  that 
°xi  =  °X2  =  ~P]  and  ex]  t  eX2*  T*11S  leads  t0  the  following  results 

for  the  linear  elastic  case: 


-exl/Pl  =  [(oi-|a2  -  L-^ L2 )  U2  +  ^  (ajj  -  ) ]/D 

Eyi / Pi  =  ^2^a2*“l  ”  *“2 ) /D 

°xl  =  _P1 

-ay1/p1  =  [L2  l£(a*  -  L^>  +  Lt  ^  («1  -  L2)]/D 


(2.53) 

(2.54) 

(2.55) 

(2.56) 
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(2.57) 


-(axl  -  ayi  )/Pi  =  ^2^ V1  +  L2l/1^V2^D 

=  -2y1cxl(l  +  ^  eff)/P] 


where 


D  =  “2^2 ^“l  "  L^)  +  “i (“2  "  L2) 
v  =  Poisson's  ratio 

"  ot1^2)/[^2^otl“2  “  +  ^Ma2  “  ^2^  (2.58) 

The  analogous  quantities  for  component,  2  are  obtained  by  interchanging  the 
indices  "1"  and  "2",  except  on  p-| .  Under  the  assumed  conditions  the  maxi¬ 
mum  resolved  shear  stress,  (a  -  a  )/ 2,  is  the  same  for  both  components. 

x  y 

The  effects  of  changing  the  boundary  condition  are  best  displayed  by 
comparing  values  of  from  Eqs.  (2.47)  and  (2.58): 


X  = 

0  .25 

.50 

.75 

1.00 

g-p-f  (Eq.  2.47) 

.0344  .0341 

.0334 

.0316 

0 

?2  eff  (^*  2.58) 

-.5932  -.590 

-.583 

-.563 

0 

yl 

-  .0116  Mb,  v>2  = 

.452  Mb 

1 

-  polyethylene,  2  - 

quartz 
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PART  III 


AN  ITERATIVE  PROCEDURE  TO  ESTIMATE  THE  VALUES  OF  ELASTIC 
CONSTANTS  OF  A  CUBIC  SOLID  AT  HIGH  PRESSURES  FROM  SOUND 
WAVE  VELOCITY  MEASUREMENTS 
Dattatraya  P.  Dandekar 


3.1  INTRODUCTION 


An  investigator  attempting  to  determine  the  variation  of  elastic  con¬ 
stants  of  solids  with  pressure  by  ultrasonic  measurements  on  new  (or  even 
well  known)  materials  may  find  that  the  needed  compressibility  measurements 
are  either  unavailable  or  if  available  are  unreliable.  Cook's  method 
enables  one  to  obtain  an  estimate  of  the  values  of  the  elastic  constants  of 

a  solid  at  high  pressure  without  a  priori  knowledge  of  the  compressibility 
3  1 

of  the  substance.  ’  In  developing  the  estimating  procedure  Cook  assumed 

that  the  parameter  a(P)  (c.f.  General  Notation  and  Analysis,  Eq.  (3.5)), 

remained  constant  with  pressure.  The  value  of  A(P)  at  any  pressure  P  is 

3  2 

given  by  its  magnitude  at  one  atmosphere.  Ruoff  '  extended  the  results 
of  Cook  in  the  case  of  cubic  solids  by  presenting  an  estimating  procedure 
which  permitted  the  parameter  A(P)  to  vary  with  pressure.  This  was  done 
by  expressing  a(P)  in  a  power  series  expansion  given  by  (3.1): 


A(P) 


a(P=1 )  +  P 


fSA(P)^ 

■+  *— 

’s2a(P)' 

8P 

V,  ^ 

m  2 

U2  J 

P+1 


+  .  .  . 


(3.1) 


where  the  quantities  on  the  right-hand  side  of  (3.1)  are  evaluated  at  one 
atmosphere. 


Even  so  the  lack  of  relevant  data  in  the  case  of  most  materials  limits 
one  to  the  first  derivative  of  a(P).  This  is  easily  seen  by  differentiating 
a ( P )  with  respect  to  pressure  P.  The  present  work  develops  an  iterative 
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procedure  to  estimate  the  values  of  the  elastic  constants  of  cubic  solids 
at  high  pressure  which  differs  from  the  one  developed  by  Ruoff  with  respect 
to  the  assumptions  regarding  (i)  the  pressure  derivative  of  the  thermal 
volume  expansion  coefficient  at  a  temperature  T,  (ii)  the  temperature 
derivative  of  the  volume  thermal  expansion  coefficient  at  a  pressure  P,  and 
(iii)  the  estimation  procedure  for  a(P).  It  is  shown  here  that  no  assump¬ 
tions  regarding  (i)  and  (ii)  are  necessary  in  order  to  estimate  the  elastic 
constants  of  cubic  solids  at  higher  pressures  provided  the  ultrasonic  measure¬ 
ments  are  made  as  a  function  of  pressure  at  more  than  two  temperatures.  This 
enables  one  to  compute  a  more  realistic  estimate  of  elastic  constants  of  cubic 
solids  as  a  function  of  pressure. 


The  size,  density  and  elastic  constants  of  a  material  specimen  change 
with  the  application  of  pressure.  Concomitant  changes  are  observed  in  the 
value  of  the  resonant  or  null  frequencies  of  a  standing  wave  and  also  in 
the  measurement  of  travel-time  for  a  pulse  between  flat  parallel  faces  of  the 
specimen.  The  analysis  presented  in  this  paper  refers  to  frequency  measure¬ 
ments  but  is  equally  valid  for  the  travel -time  measurements  of  an  elastic 
wave  propagated  in  a  medium. 

3.2'  GENERAL  NOTATION  AND  ANALYSIS 

By  a  solid  we  always  refer  to  a  cubic  solid.  Even  though  the  quan¬ 
tities  dealt  with  here  refer  to  a  pressure  P  and  a  temperature  T,  for 
simplicity  the  relevant  suffix  for  the  temperature  is  dropped  from  the 
general  notation. 


p(P): 
e(P): 
Cp(P) : 

BS(P) : 
BT(P): 


the  density  of  the  material  at  pressure  P 

volume-expansion  coefficient  of  the  material  at  pressure  P 

specific  heat  at  constant-pressure  of  the  material  at 
pressure  P 

adiabatic  bulk  modulus  of  the  material  at  pressure  P 
isothermal  bulk  modulus  of  the  material  at  pressure  P 
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XT(P):  isothermal  compress ibil ity  of  the  material  at  pressure  P 

L(0,P):  the  thickness  of  the  specimen  used  in  the  measurement  of 

the  Jth  velocity  mode  at  pressure  P 

L(J,P1 ) 

MP)  -  |_ ( j ^ p)  :  P-j  <  Pj  P-j  =  1  =  1  atm  or  1  bar,  only  in  the  case  of 

cubic  material 

V(J,P):  the  Jth  velocity  mode  in  the  material  at  pressure  P 

t(J,P):  the  travel-time  for  the  Jth  velocity  mode  at  pressure  P 

F(I,J,P):  the  Ith  null  frequency  observed  for  the  Jth  velocity  mode 

in  the  material  at  pressure  P 

N(I,J,P):  the  number  of  1/2  wavelengths  in  the  specimen  corresponding 

to  F(I,J,P) 

t(I,J,P):  the  travel-time  in  the  specimen  corresponding  to  F(I,J,P) 

IMP(J,P):  mechanical  impedance  of  quartz  transducer  for  Jth  velocity 

mode  at  pressure  P 

K(I,J,P):  IMP(J ,P)/mechanical  impedance  of  the  material  corresponding 

to  t(I,J,P) 

V(1,P):  longitudinal  velocity  in  the  (100)  direction  at  pressure  P 

V(2,P):  shear  velocity  in  the  (100)  direction  at  pressure  P 

V(3,P):  longitudinal  velocity  in  the  (110)  direction  at  pressure  P 

We  need  only  know  any  three  independent  velocity  modes  in  order  to  ob¬ 
tain  the  three  elastic  constants  of  a  solid.  In  this  paper  the  resonant 
frequencies  measured  as  a  function  of  pressure  for  the  longitudinal  modes  of 

propagation  in  the  (100)  and  (110)  directions  and  the  shear  mode  of  propaga- 

3  3 

tion  in  the  (100)  direction  have  been  used.  ’ 

We  also  assume  the  following: 

(i)  The  temperature  dependence  of  the  volume,  or  the  linear  expansion  coef¬ 
ficient  at  a  temperature  T  and  one  atmosphere  is  known; 

(ii)  the  specific  heat  at  temperature  T  and  one  atmosphere  is  known;  and 


(iii) 


3B(P) 

1  ■ 

'ae(P1)' 

3T  J 

P 

3T  J 

where  P  P-j  >  holds 


3.4 


1 


Then  the  procedure  outlined  below  can  be  used  to  estimate  the  elastic  con¬ 
stants  of  solids  at  higher  pressures,  without  reference  to  a  priori  knowledge 
of  the  compressibility  of  the  substance. 

O 

The  relation  between  the  adiabatic  bulk  modulus  and  V  (J,P),  (J  =  1,3), 
in  a  cubic  solid  may  be  written  as 

BS(P)  =  (1/3)  x  P(p){4  x  V2(3,P)  -  4  x  V2(2,P)  -  V2(1,P)}  (3.2) 

Relation  (3.1),  expressed  in  terms  of  L(J,P-j),  x(J,P),  x(P)  and  p(P-|),  is 
given  as  relation  (3.3): 

s  4  X  L2(3,P,)  4  x  L2(2 ,P, )  L2(l ,P, ) 

Bb(P)  =  (l/3)p(P, )  x  X(P)  x  {- - = - 1 - * - - = - '—}  (3.3) 

T  (3,P)  t2(2,P)  t2(1,P) 

where  p(P)  =  x3(P)  x  p(P1). 

By  the  definition  of  isothermal  bulk  modulus  we  obtain 


B 1 (P)  =  -Vol .  (P) 


where  Vol.  (P)  =  l/p(P).  If 


aP  ] 

=  p(P) 
T 

3P  1 

-  A(P) 

3P  ] 

(9  Vol.  (P)J 

L'ScTpTJ 

T  3 

[3X(P)J 

(3.4) 


a(P) 


g2(P)  x  BS(P)  x  T 
p(P)  x  Cp(P) 


where  temperature  T  is  in  degrees  Kelvin,  then 

rJ,n  _  BS(P) 

B  (p)  ~  rr+Tpyj' 


(3.5) 


(3.6) 


3  5 

Using  Williams  and  Lamb's  *  method  of  ultrasonic  wave  velocity  measurements 

O  C 

as  modified  by  Colvin,  transit  time  for  the  various  wave  propagations  is 
obtained  from  the  following  relations: 
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(3.7) 


N(I,J,P)  =  Integer^*  -  0.5  -  K(I,J,P)} 


N(I,J,P)  +  0.5  K(I,J,P)r  1  1  i 

2  x  F(I,J,P)  "  2  <F(R,J,P)  '  FTOTPl1 


(3.8) 


In  the  above  expressions  K(I,J,P)  may  be  written  as 


K(I,J,P) 


,  IMP(J?P)  T  =  -IMP(J.P)  X  T(J,P) 
p(P)  x  V(J>)  p(Pi)  x  x  l(J,P1) 


(3.9) 


where  IMP(0,P)  is  the  mechanical  impedance  of  the  transducer  for  the  Jth 
velocity  mode  at  pressure  P. 


It  is  evident  from  relation  (3.8)  that  if  the  measurements  are  made 
near  F(R,J,P)  any  error  in  the  estimation  of  t(I,0,P)  due  to  inaccurate 
knowledge  of  K(I,J,P)  becomes  negligible. 


By  integrating  relation  (3.4)  we  obtain 


X(P) 


=  x(P,)  x  exp. 


f  P  -  P, 


3  x  BT(P) 


(3.10) 


Two  things  should  be  noted  regarding  the  derivation  of  (3.10)  from  (3.4): 

(i)  in  the  definition  of  isothermal  bulk  modulus  at  a  pressure  P,  one  could 
obtain  its  value  by  either  decreasing  or  increasing  the  pressure  slightly; 
and  (ii)  when  integrating  (3.4)  it  must  be  remembered  that  it  is  implied 
in  the  definition  of  B^(P)  that  it  remains  constant  over  the  range  of  inte¬ 
gration  P-|  to  P.  In  expression  (3.10)  it  is  implied  that  the  isothermal 
bulk  modulus  of  a  substance  at  pressure  P  has  been  obtained  by  decreasing  the 
pressure  from  P  to  P-j .  The  expression  for  A(P)  as  derived  above  differs  from 
that  obtained  by  following  either  Cook's  or  Ruoff's  procedures.  The  expres¬ 
sion  for  a(P)  that  will  be  obtained  by  following  Cook's  or  Ruoff's  procedure 
may  be  given  by 
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A(P)  =  1  +  - '-2 -  f,  [1  +  A(P)] 

p(D  L2(1)  1 


t2(3,P)  t2(2,P)  t2(1,P) 


(3.11) 


where  L ( J ,  1 )  =  L(l),  and  [1  +  a( P) ]  is  a  constant  in  Cook's  method  and 
equals  [1  +  a(1)]  but  is  a  variable  in  Ruoff's  method. 

This  expression  (3.11)  for  x(P)  is  arrived  at  by  expressing  BT(P)  in 
terms  of  p(l),  L ( 1 ) ,  t(J,P),  and  a(P)  with  the  help  of  Eqs.  (3.3)  and  (3.6), 
and  integrating  (3.4).  Hence,  \(P)  in  (3.11)  can  be  determined  if  the  value 
of  a(P)  can  be  estimated.  Ruoff  *  estimates  the  value  of  a(P)  from  rela¬ 
tion  (3.1)  by  means  of  thermodynamic  relations.  For  example,  to  evaluate 

C 

( 3A( P)/3P)j ,  we  would  rewrite  B  (P)  as 


BS(P)  =  p(P)  x 


t2(P) 


(3.12) 


where  L2(P)/t2(P)  is  an  abbreviation  for  the  coefficient  of  p(P)  in  Eq.  (3.2). 
Then  a(P)  in  (3.5)  may  be  written  as 


a  ( p )  =  »-T.  jl.&pJ 

T  (P)  X  Cp(P) 

and  the  logarithmic  derivative  of  a(P)  yields 


(3.13) 


1  f3A(P)l  _  2  f3g(P)1  .  2  faL(P)' 

wn  aP  JT  "  WJ{  9P  JT  lTpTI  9P 


2  fax(P)|  1  f9Cp(p^ 
7[py[  3P  JT  ■  C^PTI,  3P 


(3.14) 


From  thermodynamics,  we  know  that 


3B(P) 

3P 


(3.15) 
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and 


faCp(P)] 


3g(P) 


9T 


+  3^(P)} 


Hence,  in  the  limit  as  P  1,  the  expression  (3.14)  reduces  to 


(3.16) 


1 

9A(P) 

2 

bxV>] 

aOT 

9P  J 

T,P+1 

_  3HT 

9T  J 

2  T 

3  x 


(1) 


am 

9P 


'  T ,  P-+1 


Til 


c^pT 


<(■ 


3g(P)l 


9T 


P+1 


+  32(D}. 


(3.17) 


Thus  the  magnitude  of  the  first  derivative  of  a(P)  in  the  limit  as  P 1  may 
be  determined  if  (9x^(P)/8T)p_i ,  (9p(P)/9T)p_-| ,  p(l),  and  Cp ( 1 )  are  known  a 
and  the  value  of  a(P)  may  be  approximated  at  a  pressure  P  by 


a(P) 


a(1)  +  P 


fMiP) 

9P 


'P-+1  ,T 


(3.18) 


Similarly  the  higher  derivatives  of  a(P)  may  be  evaluated  if  the  relevant 
thermodynamic  data  are  available. 


The  expression  for  x(P)  in  the  new  method,  i.e.  relation  (3.10)  des- 
described  in  this  paper,  is  seen  to  differ  from  the  earlier  two  works  for 
two  reasons.  In  their  works,  (i)  x(P)  is  defined  as  L(J ,1 )/L(J ,P) ,  and 
(ii)  a(P)  is  estimated  by  a  different  procedure. 


The  quantities  measured  or  known  are  p(l),  L(J,1),.  F(I,J,P)  or  t(J,P), 
P  and  T.  For  quartz  transducers,  IMP(J,P)  and  F(R,J,P)  can  be  obtained 
safely  to  4000  bars  and  from  room  temperature  to  90° K  from  the  work  of 
McSkimin  and  Andreatch.  '  This  information  is  not  required  if  the  ultra¬ 
sonic  measurements  are  of  the  travel-times.  g(P)  is  usually  known  only  as 
a  function  of  temperature  at  one  atmosphere.  However,  the  variation  in  the 
elastic  constants  with  temperature  at  pressure  P  provides  one  with  the  tem¬ 
perature  derivative  of  the  isothermal  compressibility.  And  from  relation 
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(3.15)  one  may  obtain  b(P)  at  temperature  T  if  $(P)  is  known  at  one  atmo¬ 
sphere  and  temperature  T.  In  a  normal  substance  where  P^  £  P, 


kV)] 

'3X(P1)1 

3T 

\  J 

< 

P 

3T  J 

holds.  So,  to  assume  that 


f3p(P)l 

kV,)] 

l  &  J 

T 

aT 

(3.19) 


(3.20) 


ensures  that  the  value  of  &(P)  obtained  from  (3.5)  is  underestimated.  If 
ultrasonic  measurements  are  made  as  a  function  of  pressure  at  more  than  one 
temperature,  a  better  estimate  of  (3x^(P)/3"Op  may  be  obtained  by  simply 
incorporating  (3x^(P)/3T)p  as  an  additional  parameter  to  be  iterated  according 
to  the  scheme  presented  in  Figure  3.2.  Where  such  information  is  unavail¬ 
able  (3.19)  or  (3.20)  may  be  used.  Similarly  the  computation  of  Cp(P)  from 
relation  (3.16)  by  assuming 


3B(P) 

'3B(P-|)' 

3T  j 

T 

3T 

V.  / 

(3.21) 


implies  that  the  resulting  values  of  Cp(P)  from  relation  (3.16)  will  also 
be  underestimated.  However,  the  resultant  error  in  the  estimated  value  of 
A(P)  due  to  the  intrinsic  underestimation  of  e(P)  and  Cp(P)  is  likely  to  be 
small,  up  to  3-4  kilobars  for  most  materials.  Thus  everything  in  expression 

C 

(3.5)  except  B  (P)  and  p(P)  is  either  known  or  may  be  approximated  with 
reasonable  accuracy. 


The  iterative  procedure  described  below  is  that  presented  in  Figure  3.1, 
because  we  feel  that  the  understanding  of  the  procedure  given  in  Figure  3.2 
will  be  facilitated  by  an  understanding  of  the  simpler  procedure.  Thus  the 
iterative  procedure  described  assumes  that  relations  (3.20)  and  (3.21)  hold. 


At  P  =  1  atmosphere,  all  the  quantities  involved  are  known;  no  itera¬ 
tion  is  required  to  estimate  the  required  elastic  constants  of  solids. 
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FIGURE  3.1 

A  Flow  Chart  of  the  Iterative 
Procedure  to  Estimate  the  Varia¬ 
tion  in  the  Elastic  Constant  of 
a  Cubic  Solid  with  Pressure  When 
the  Elastic  Wave  Velocities  are 
Obtained  from  the  Measurement  of 
the  Resonant  Frequencies  of  a 
Standing  Wave  as  a  Function  of 
Pressure  at  a  Fixed  Temperature. 


If  X'(P)=X(P) 

jk _ 

All  other 
parameters 
at  pressure  P 


43 


Increase  Pressure  (P) 

^  < - - - 

Temperature  (T)  ^ _ 

^  next  temperature 

Set  A(P,T)=A(Preceding  Pressure,  T) 


Set  A(P,T)-A'(P,T) 


B(P,T) 


Set  x(P,T)=x(Preceding  Pressure,  T) 


Set  x(P,T)=x' (P,T) 


Jth  velocity  mode  at  ^ - 7 — : — rr - T~ 

Pressure  P  and  Temperature  T  next  ve^oc1ty  moc*e 


KJ.P.T) 

V(J,P,T) 


When  all  velocity  modes 
are  calculated  at  Pressure  P 
and  Temperature  T 


BS(P,T) 

aCp(P.T)1 

aP  ' 
a(P,T) 
BT(P,T) 
X 1 (P,T) 


If  X' (P,T)/x(Preceding  Iteration,  T) 


If  x' (P,T)=x(P,T) 


When  all  B^(P,T)  are 
calculated  at  Pressure  P 

, _ ±— 

A'(P.T)  I 


If  A' (P,T)^A(Preceding  Iteration,  T) 


All  other  parameters  at 


FIGURE  3*2  If  A'(P,T)=A(P,T) 

A  Flow  Chart  of  the  Iterative 

Procedure  to  Estimate  the  - - - — 

Variation  in  the  Elastic  Con-  All  other  parameters  at 

stants  of  a  Cubic  Solid  with  Pressure  P  and  all  Temperatures 

Pressure  when  the  Travel -time  L _ _ _ _  _ K  _ 

Measurements  are  made  as  a 

Function  of  Pressure  at  more  ^ - - — — - 

than  Two  Temperatures.  A(P,T) 

=  [3x(P»T)/3T]p.  End 
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At  the  next  higher  pressure  all  the  fundamental  quantities  in  the 
relations  (3.3),  (3.5),  (3.6),  (3.7),  (3.8),  (3.9),  and  (3.10),  except  A (P) 
and  K(I,J,P)  are  known.  The  procedure  developed  here  involves  a  two  stage 
iteration,  one  at  the  level  of  pressure  and  the  other  on  the  Ith  null  fre¬ 
quency  of  the  Jth  mode.  We  set  A (P)  =  A(preceding  pressure)  and  K(I,J,P)  = 

1  and  estimate  N(I,J,P)  and  t(I,J,P)  and  K(I,J,P).  If  the  value  of  K(I,J,P) 
thus  obtained  agrees  with  the  previously  assigned  value  we  compute  N(I,J,P) 
for  the  (I  +  1 ) th  frequency.  If  this  value  of  K(I,J,P)  does  not  agree  with 
the  previously  assigned  value  these  values  of  N(I,J,P)  and  t(I,J,P)  are  cor¬ 
rected  by  setting  K(I,J,P)  equal  to  the  value  obtained  last,  and  iterating 
all  over  again.  This  is  repeated  till  two  consecutive  estimates  of  K(I,J,P) 
are  the  same.  A  similar  computation  is  performed  for  all  the  velocity  modes. 
By  interpolation,  from  these  t(I,J,P)'s  one  obtains  values  corresponding  to 
F(R,J,P),  each  of  which  is  called  x(J,P).  These  t(J,P)'s  in  turn  are  used 
to  obtain  V(J,P)  which  together  with  p(P)  yield  an  estimate  of  B^(P),  a(P), 
B^(P),  and  finally  a(P).  If  the  value  of  \(P)  thus  obtained  agrees  with  the 
previously  assigned  value,  the  estimates  of  B  (P)  and  B  (P)  are  correct.  If 
this  value  of  a(P)  does  not  agree  with  the  previously  assigned  value  these 
B  (P)  and  B  (P)  are  corrected  by  means  of  setting  a(P)  equal  to  the  value  of 
A(P)  obtained  last  and  iterating  all  over  again.  This  is  repeated  till  two 
consecutive  estimates  of  \(P)  are  the  same.  Once  this  is  known  all  other 
elastic  constant  parameters  may  be  obtained.  This  iterative  procedure  is 
sketched  diagramatical ly  in  Figure  3.1. 

Table  3.1  displays  the  estimates  of  the  pressure  derivatives  of  the 
adiabatic  and  isothermal  bulk  moduli  of  NaCl  and  KC1  at  295°  and  195°K  ob¬ 
tained  from  the  above  mentioned  iterative  procedure.  The  required  travel 
time  data  as  a  function  of  pressure  for  this  computation  were  reconstructed 

from  the  pressure  derivatives  of  the  travel -time  for  the  various  elastic  wave 

3  ft 

velocities  given  in  the  paper  of  Bartels  and  Schuele.  All  other  ancillary 
data  used  were  also  taken  from  reference  3.8.  It  may  be  seen  that  the  esti¬ 
mates  of  the  pressure  derivatives  of  the  bulk  moduli  of  NaCl  and  KC1  obtained 
in  the  present  work  for  pressures  ranging  up  to  1.7  kb  differ  slightly  from 
those  obtained  by  Bartels  and  Schuele.  However,  such  differences  may  become 
significant  at  higher  pressures.  It  should  be  noted  further  that  the  iterative 
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TABLE  3.1 


The  Pressure  Derivative  of  the  Adiabatic  and  Isothermal  Bulk  Moduli  of 
NaCl  and  KC1  as  Obtained  by  Bartels  and  Schuele  (B  and  S),  and  as 
Obtained  in  the  Present  Work  (D)  from  the  Data  of  Bartels  and  Schuele. 


295°K 

Adiabatic 

5.27 

5.33 

Isothermal 

5.35 

5.38 

1 95°K 

Adiabatic 

5.13 

5.18 

Isothermal 

KC1 

5.20 

5.23 

295°K 

Adiabatic 

5.34 

5.36 

Isothermal 

5.41 

5.44 

195°K 

Adiabatic 

5.34 

5.36 

Isothermal 

5.41 

5.43 
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procedure  outlined  in  the  present  work  may  be  easily  applied  to  determine 
the  variation  in  the  elastic  constants  of  an  isotropic  solid. 

Figure  3.2  is  the  schematic  representation  of  the  iteration  procedure 
when  the  travel -time  measurements  are  made  as  a  function  of  pressure  at  more 
than  two  temperatures. 

We  are  in  the  process  of  developing  a  variant  of  this  iterative  pro¬ 
cedure  designed  to  estimate  the  elastic  constants  of  a  non-cubic  solid  as 
a  function  of  pressure. 


47 


REFERENCES 


3.1  R.  K.  Cook,  J.  Acoust.  Soc.  Am.  29_,  445  (1957). 

3.2  A.  L.  Ruoff.,  Intra-Laboratory  Report  No.  8,  Thurston  High  Pressure 
Laboratory,  Cornell  University,  Ithaca,  New  York  (1968). 

3.3  Change  in  the  choice  of  independent  velocity  modes  merely  changes 
the  expression  for  the  adiabatic  bulk  modulus  in  expression  (3.2). 

3.4  Assumption  (iii)  need  not  be  made  if  ultrasonic  measurements  are 
made  as  a  function  of  pressure  at  more  than  two  temperatures. 

3.5  J.  Williams  and  J.  Lamb,  J.  Acoust.  Soc.  Am.  30^,  308  (1958). 

3.6  A.  D.  Colvin,  M.Sc.  Thesis,  Rensselaer  Polytechnic  Institute,  Troy, 

New  York  (1959). 

3.7  H.  J.  McSkimin  and  P.  Andreatch,  J.  Acoust.  Soc.  Am.  3^,  609  (1962). 

3.8  R.  A.  Bartels  and  D.  E.  Schuele,  J.  Phys.  Chem.  Solids  26_,  534  (1965). 


48 


PART  IV 


AN  ITERATIVE  SCHEME  TO  ESTIMATE  THE  VALUES  OF  ELASTIC  CONSTANTS 
OF  A  SOLID  AT  HIGH  PRESSURES  FROM  SOUND  WAVE  VELOCITY  MEASUREMENTS 

Dattatraya  P.  Dandekar 


4.1  INTRODUCTION 

In  Section  III  we  presented  an  iterative  procedure  to  estimate  the 
values  of  elastic  constants  of  cubic  and  isotropic  solids  at  high  pressures 
from  sound  wave  velocity  measurements  under  the  assumption  that  the  concomi¬ 
tant  compressibility  measurements  are  unknown.  This  procedure  resembled 

4  1  4  2 

procedures  developed  by  Cook  '  and  Ruoff  '  with  regard  to  the  use  of  the 
principle  of  self-consistent  integration  but  differed  from  both  with  regard 
to  the  assumption  of  the  variation  in  a(P,T),  defined  as  &  (P,T)  x  T/x  ( P ,T) 
x  p ( P ,T)  x  Cp(P,T),  with  pressure  and  temperature  and  also  with  regard  to 
the  method  of  estimating  a(P,T).  This  work  presents  a  similar  procedure 
which  permits  estimation  of  the  values  of  elastic  constants  of  a  noncubic 
crystalline  solid  at  high  pressures  from  measurements  of  sound  wave  velo¬ 
cities.  The  iterative  scheme  presented  here  is  a  general  one.  When  the 
ultrasonic  measurements  are  made  as  a  function  of  pressure  at  more  than  two 
temperatures  no  restrictive  assumptions  are  required.  However,  when  such 
measurements  are  made  as  a  function  of  pressure  at  only  one  temperature,  two 
additional  assumptions  become  necessary  before  the  iterative  procedure  can 
be  applied. 

These  two  assumptions  are: 

(i)  The  temperature  derivatives  of  the  linear  compressibilities  of  a  solid 
are  independent  of  pressure. 

(ii)  The  temperature  derivatives  of  the  linear  thermal  expansion  coeffi¬ 
cients  are  independent  of  pressure. 

If  measurements  are  made  at  two  temperatures,  assumption  (ii)  is  required. 
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For  simplicity  of  presentation  the  analysis  presented  refers  to  mea¬ 
surements  of  the  transit  time  of  an  elastic  wave  propagated  in  an  anisotropic 
medium  rather  than  to  measurements  of  resonant  frequency.  Since  application 
of  pressure  changes  the  density,  size  and  elastic  properties  of  a  solid, 
which  is  reflected  by  a  concomitant  change  in  the  transit  time  of  an  elas¬ 
tic  wave  propagated  in  the  solid,  the  change  in  transit  time  with  pressure 
contains  all  the  desired  information  about  the  changes  brought  about  in  a 
solid  due  to  the  application  of  pressure. 


4.2  CONVENTIONS  AND  NOTATIONS 


The  elastic  constants  of  a  solid  may  be  represented  by  reference  to 

any  Cartesian  coordinate  system  (x,y,z).  However,  it  is  convenient  to  refer 

them  to  a  system  of  coordinates  having  some  relation  to  the  crystal  axes  < 

(a,b,c).  In  this  work,  whenever  applicable  the  relationship  between  the 

Cartesian  coordinate  system  (x,y,z)  and  the  crystal  axes  (a,b,c)  is  the  one 

4  3 

adopted  by  the  IRE  Standards  Committee.  '  The  number  of  elastic  constants 
(Ci j kt )  necessary  to  characterize  the  elastic  property  of  a  solid  depends 
upon  the  crystal  class  to  which  the  solid  belongs.  The  subscripts  of  these 
constants  were  contracted  to  (C  )  by  following  the  usual  convention  of 
writing  the  subscript  ij,  kt  =  11 ,  22,  33,  23,  12,  12  by  p,  q  =  1,  2,  3, 
4,  5,  6  so  as  to  represent  these  Cp^  by  a  6x6  matrix,  denoted  by  [Cp  ].  The 
corresponding  6x6  matrix  of  the  elastic  compliances  Sp^  may  be  obtained  from 
Cpq  by  using  the  matrix  relation  between  them,  namely 


[Cpq]x[Spq]  -  [I]  (4.1) 

where  [  I  ]  is  a  6x6  unit  matrix.  The  C  matrices  for  different  crystal 

pq 

classes  may  be  found  in  reference  4.4. 


4.2.1  Notation 


P: 

T: 


Pressure 

Temperature 
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p(P,T):  Density  of  a  solid  at  pressure  P  and  temperature  T 

x(£,m,n,P,T) :  Linear  compressibility  of  the  solid  in  the  direction  whose 

direction  cosines  are  given  by  a,  m,  and  n  at  P  and  T 
3(«,,m,n ,P,T) :  Linear  thermal  expansion  coefficient  of  the  solid  in  the 

direction  whose  direction  cosines  are  given  by  a,  m,  and 
n  at  P  and  T 

x(PjT):  Volume  compressibility  of  the  solid  at  P  and  T 

e(P,T):  Volume  thermal  expansion  of  the  solid  at  P  and  T 

L(t,m,n,J,P,T):  Width  of  the  specimen  used  to  measure  the  Jth  velocity  mode 
in  the  solid  in  the  direction  such  that  i,  m,  n  are  the 
direction  cosines  at  P  and  T 

-rU,m,n,J,P,T):  Transit  time  of  the  Jth  wave  velocity  mode  corresponding  to 
LU,m,n,J,P,T) 

VU,m,n,J,P,T) :  The  Jth  velocity  in  the  solid  in  the  direction  whose  direc¬ 
tion  cosines  are  a,  m,  and  n  at  P  and  T 
LU,m,n,J,P,,T) 

x(£,m,n,J,P,T):  -  n^-  y-j ^  -  where  P^  <  P,  and  P-|  equal  to  unity  indi 

cates  one  atmosphere  pressure 

Cp(P,T):  Specific  heat  of  the  solid  at  constant  P  and  T 

The  superscripts  T  and  S  attached  to  a  quantity  indicate  its  isothermal  and 
adiabatic  value  respectively. 

4.3  PRELIMINARIES 

The  various  relations  used  in  the  iterative  procedure  to  be  described 
originate  either  in  the  theory  of  elastic  wave  propagation  in  a  solid  or  in 
thermodynamic  theory.  These  general  relations  are  presented  with  brief  intro¬ 
ductory  remarks  drawn  from  the  two  theories  and  appended  only  to  clarify  the 
material  of  this  paper. 

4.3.1  Relations  Obtained  from  the  Theory  of  Elastic  Wave  Propagation 
The  elastic  constants  of  a  solid  are  determined  by  measuring 


the  velocities  with  which  elastic  waves  are  propagated  along  several  direc¬ 
tions  in  the  solid.  The  number  of  velocity  measurements  needed  to  understand 
the  elastic  property  of  the  solid  depends  on  the  crystallographic  class  to 
which  it  belongs.  Christoffel's  equations  (e.g.  Eq.  (4.2)),  which  are 
applicable  to  any  crystalline  system,  give  the  relationship  between  measured 
velocities  and  elastic  constants.  In  general,  for  a  plane  wave  propagated 
in  a  crystal  having  direction  cosines  i,  m,  n,  the  three  possible  wave  velo¬ 
cities  V  may  be  found  in  terms  of  the  elastic  constants  C  from  the  roots  of 

pq 

Christoffel's  equations. 


A11  ‘  py2  A1 2 

A12  AA22  -  pV2 

A13  A23 


A23 

A23 

A33  "  pV 


(4.2) 


It  follows  that  p V  is  related  to  C.. in  a  manner  which  is 
determined  by  the  direction  in  which  a  wave  is  propagated.  Three  different 


velocities  of  propagation  imply  that  the  three  displacement  vectors  associ¬ 
ated  with  these  velocities  are  mutually  perpendicular  and  hence  independent. 
Usually  the  three  waves  are  mixed;  one  is  predominantly  longitudinal  and  the 


other  two  are  predominantly  shear.  Pure  waves  may  be  propagated  only  in  a 
few  special  crystallographic  directions.  The  location  of  the  pure  mode 

directions  in  crystals  of  various  symmetries  have  been  investigated  exhaus- 

4  5  4  6 

tively  by  Borgnis  *  and  Brugger.  ' 


4.3.2  Thermodynamic  Relations 

These  relations  serve  two  purposes: 
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(i)  To  evaluate  the  pressure  derivatives  of  the  specific  heat  and  the 
linear  thermal  expansions  in  the  three  principal  directions  of  the 
solid  at  pressure  P  and  temperature  T,  and 

(ii)  To  convert  the  adiabatic  quantities  into  their  isothermal  counterparts 
at  pressure  P  and  temperature  T. 

The  adiabatic  and  isothermal  elastic  compliances  are  related  by 

sijkt(P'T)  '  s1jkt(p’T)  -  -sij(P,T)  6kt(PJ)  '  T/cp(p’T)  x  p(p,T)  (4’4) 

When  expressed  in  our  notation  the  following  relation  between  x  U,m,n,P,T) 
and  xTU>ro>n>pJ)  results. 

X  U,m,n,P,T)  =  x  ( £»m,n »P»T)  +  *x^Cp"("P,T) 

{3(1 ,0,0,P,T)jt2  +  6(0,1 ,0,P,T)m2  +  e(0,0,l  ,P,T)n2}  (4.5) 

where 

6(P,T)  ee  6(1,0, 0, P,T)  +  6(0, l.O.P.T)  +  6(0,0,1,PJ)  (4.6) 

and 

xU,m,n,P,T)  =  (S^-|  +  S12  +  S13)  2  +  (S12  +  S22  +  S23)m2 

+  (S13  +  S23  +  S33)n  (4.7) 

I 

Equation  (4.5)  may  also  be  written  as 

xT(£,m,n,P,T)  =  xSU>m,n,P,T)  +  A(£,m,n,P,T)  (4.8) 

where 

A(^m,n,P,T)  =  p(p^>IJCp<P,T)  <e(1  ’0-°-P-T)  * 

+  6(0 ,1 ,0 »P ,T)m  t  3 ( 0 ,0 ,1 , P ,T ) n  } 
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From  the  definition  of  isothermal  linear  compressibility,  we  have 


v^((  m  p  p  t]  =  _  1 _ r  ExX] 

X  ^,m,n,K,u  L()i,m,n,P,T)  I  3P 


1 _  [3xU,m,n,P,Tl 

x(£,m,n,P,T)  I  3P 


(4.9) 


Since  the  subscript  J  has  no  significance  in  the  above  relation,  it  has  been 
omitted  from  L(t,m,n,J,P,T) .  By  integrating  Eq.  (4.9)  with  respect  to  pres¬ 
sure,  we  obtain 

x(£,m,n,P,T)  =  xU,m,n,P-|  ,T)  x  exp{(P-P])  xTU,m,n,P,T)}  (4.10) 

The  above  integration  is  performed  by  assuming  that  the  compressibility  is 
determined  by  reducing  the  pressure  from  P  to  P,  and  by  taking  into  account 
the  definition  (4.8)  which  implies  x  U,m,n,P,T)  remains  constant  in  the 
range  of  integration  P  to  P-| .  Again  the  temperature  derivative  of  linear 
compressibility  is  related  to  the  pressure  derivative  of  the  linear  thermal 
expansion  of  a  material  by 


3x  (£,m,n,P,T! 
3T 


3gU,m,n,P,T] 

3P 


(4.11) 


and  the  pressure  derivative  of  specific  heat  may  be  written  as 


3Cp(P,T) 


JTpTTT  {b2(p.t)  *  1 


(4.12) 


The  use  of  the  above  set  of  relations  enables  one  to  estimate  the  values  of 
the  elastic  constants  of  a  solid  at  high  pressure  without  a  priori  knowledge 
of  the  compressibility  of  the  substance. 


4.4  GENERAL  ITERATIVE  SCHEME 


In  general  the  iterative  scheme  proposed  here  attempts  to  obtain  self- 
consistent  estimates  of  x(£,m,n,P,T)  in  the  three  principal  directions,  i.e. 
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X(1,0,0,P,T),  X(0,1,0,P,T)  and  x(0,0,l,P,T)  at  pressure  P  and  temperature 

T.  The  iterative  procedure  presented  below  assumes  the  following: 

(i)  The  temperature  dependence  of  the  linear  expansion  coefficients  are 
known  at  one  atmosphere. 

(ii)  The  value  of  specific  heat  is  known  as  a  function  of  temperature  at 
one  atmosphere. 

(iii)  An  adequate  number  of  velocity  measurements  are  made  to  extract 
information  about  elastic  properties  of  the  solid  as  a  function  of 
pressure  at  more  than  one  atmosphere. 


Then,  at  P  =  1  and  within  the  experimental  range  of  temperature,  all 
the  quantities  in  the  above  set  of  12  relations  are  known;  no  iteration  is 
required  to  estimate  the  values  of  the  elastic  constants  of  a  solid  as  a 
function  of  temperature  at  one  atmosphere. 


At  the  next  higher  pressure  all  the  quantities  except  xU,m,n,P,T)  in 
the  above  set  of  relations  are  known.  The  procedure  developed  here  involves 
a  three-stage  iteration  in  order  to  yield  a  self-consistent  estimate  of  not 
only  the  elastic  constants  but  also  of  pressure  and  temperature  derivatives 
of  the  linear  thermal  expansion  and  of  the  specific  heat  at  pressure  P  and 
temperature  T.  The  three  stages  of  iteration  are  done  at  the  levels  of 
pressure,  temperature  and  the  three  principal  directions  of  the  solid.  To 
initiate  the  procedure,  we  set 


'  ,  > 

T 

8x  U,m,n,P,T) 

T  1 

3X  (£,m,n,P] ,T ) 

8T  j 

P 

3T  J 

and 


3BU»m,n,P,T) 

'ae(z,m,n,P-|  J)y 

3T  j 

P 

3T  J 

(4.14) 


where  P-|  is  the  preceding  value  of  pressure.  These  enable  us  to  compute 
e(£,m,n,P,T)  and  Cp(P,T). 
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Next  we  set 


xU,m,n,P,T)  =  x(t,m,n,P^,T)  (4.15) 

With  this  assumption,  we  compute  p(P,T)  from  the  relation 


xT(PJ) 


=  xT(i >0,0, P,T)  +  xT(o,i,o,p,t)  +  xt(o,o,i,p,t) 


(4.16) 


and  L(£,m,n,P,T);  these  with  the  known  values  of  T'U,m,n,J,P,T)  enable  us 
to  estimate  Cp^(P,T)  from  Eqs.  (4.2).  The  use  of  relations  (4.1),  (4.4), 

(4.11)  and  (4.12)  together  with  assumptions  (4.13)  and  (4.14)  provides 
estimates  of  the  values  of  sL(P,T).  From  these  estimates  of  si  (P,T)  we 
obtain  x  (£,m,n ,P,T) ,  which  in  turn  by  relation  (4.10)  yields  a  new  esti¬ 
mate  of  x(t,m,n,P,T).  If  the  new  values  of  x(t,m,n,P,T)  in  the  three 
principal  directions  agree  with  their  respective  values  assumed  at  the 
beginning  of  the  iteration,  the  estimated  values  of  the  elastic  constants 
are  correct  and  consistent  with  the  assumptions  represented  by  relations 
(4.13)  and  (4.14).  If  these  values  of  x(&,m,n,P,T)  do  not  agree  with  the 
previously  assumed  values,  the  iteration  is  repeated  with  these  new  values 
of  xU,m,n,P,T)  as  starting  values,  and  all  the  quantities  are  recalculated. 

This  process  is  repeated  until  two  successive  estimates  of  x(t,m,n ,P,T)  in 
the  three  respective  principal  directions  are  equal  in  magnitude  at  the  pres¬ 
sure  P  and  a  temperature  T.  This  whole  iteration  procedure  is  carried  out  at 
the  pressure  P  and  at  all  the  temperatures  at  which  travel -time  measurements 
are  made.  When  all  the  elastic  constants  of  a  solid  are  estimated  in  the 
above  manner  at  the  pressure  P  as  a  function  of  temperature,  we  obtain  new 
estimates  of  [3xTU>rn,n,P,T)/3T]p  and  [a&U,m,n ,P,T)/3T]p.  If  these  new 
values  of  [3xTU,m,n,P,T)/3T]p  and  [3BU,m,n,P,T)/3T]p  in  the  three  principal 
directions  agree  with  their  respective  assigned  values  at  the  beginning  of 
the  iteration,  the  estimated  values  of  the  elastic  constants  at  the  pressure 
P  and  all  the  temperatures  are  correct.  If  these  new  values  of  [3xTU,m,n ,P,T)/3T]p 
and/or  [aeU,m,n,P,T)/3T]p  do  not  agree  with  their  previously  assigned  values. 
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the  iterative  procedure  is  repeated  with  these  new  values  as  the  starting 
point.  Such  a  procedure  is  repeated  until  two  consecutive  estimates  of 
[9xT(t,m,n,P,T)/9T]p  and  [9BU,m,n,P,T)/9T]p  agree  with  their  respective 
initial  values  in  all  the  three  principal  directions.  The  iterative  pro¬ 
cedure  described  above  is  represented  by  a  flow  chart  in  Figure  4.1. 

The  net  outcome  of  the  above-mentioned  procedure  is  to  yield  (i) 

values  of  the  elastic  constants  as  functions  of  pressure  and  temperature, 

(ii)  values  of  the  partial  temperature  and  pressure  derivatives  of  the 
linear  thermal  expansion  coefficient  as  functions  of  pressure  and  tempera¬ 
ture,  and  (iii)  the  pressure  derivative  of  the  specific  heat  as  a  function 
of  temperature.  If,  however,  the  transit  times  of  the  elastic  waves  in  a 
solid  are  measured  as  functions  of  pressure  at  only  one  temperature,  then 

values  of  the  elastic  constants  of  the  solid  may  be  estimated  at  high  pres¬ 

sures  (i)  by  interpreting  P-|  in  relations  (4.13)  and  (4.14)  to  signify  the 
preceding  pressure  at  which  these  quantities  are  known,  and  (ii)  by 
omitting  the  iterations  on  [9xTU,m,n,P,T)/9T]p  and  [93U,m,n,P,T)/9T]p  in 
the  above  scheme. 

Table  4.1  gives  two  sets  of  values  of  the  elastic  constants  of  cal- 
cite  with  pressure  at  298°K.  Both  are  based  on  the  data  of  reference  4.7, 
but  the  one  labelled  "D"  was  obtained  without  iteration,  using  Bridgman's 
value  for  compressibi 1 ity;  the  other  was  obtained  from  the  present  itera¬ 
tive  procedure  in  conjunction  with  the  required  thermodynamic  data.  The 
sources  of  these  thermodynamic  data  are  indicated  under  the  table.  These 

two  sets  of  estimates  of  C-.  of  calcite  are  in  good  agreement.  A  similar 

*  J 

calculation  performed  on  the  ultrasonic  data  of  rutile  is  shown  in  Table  4.2. 

We  are  in  the  process  of  analyzing  the  relevant  ultrasonic  data  on 
several  crystalline  solids. 
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FIGURE  4.1 


A  Flow  Chart  of  the  Iterative  Scheme  to  Esti-. 
mate  the  Values  of  the  Elastic  Constants  of  a 
Solid  as  a  Function  of  Pressure  and  Tempera-  I 
ture.  AU,m,n,P,T)  =  [3xTU,m,n,P,T)/3T]p 


and  BU,m,n,P,T)  =  [seU.m.n  ,P,T)/3T] 


TABLE  4.1 
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Adiabatic  Elastic  Constants  (in  units  of  10  dyn/cm  )  of  Calcite  with  Pressure  at  298°K  as  Obtained  by 
Dandekar^'7  (D)  and  as  Obtained  by  the  Iterative  Procedure  (IP)  from  the  same  Ultrasonic  Data. 


Pressure 

C11 

C33 

C44 

C 

66 

C14 

C13 

Kbar 

D 

IP 

D 

IP 

D 

IP 

D 

IP 

D 

IP 

D 

IP 

+0.05 

+0.05 

+0.02 

+0.07 

+0.06 

+0.33 

0.001 

14.626 

14.626 

8.531 

8.531 

3.405 

3.405 

4.328 

4.328 

-2.076 

-2.076 

5.076 

5.076 

2.0 

14.650 

14.651 

8.526 

8.527 

3.453 

3.453 

4.335 

4.336 

-2.078 

-2.080 

5.190 

5.179 

4.0 

14.674 

14.676 

8.522 

8.522 

3.468 

3.468 

4.342 

4.343 

-2.101 

-2.100 

5.419 

5.402 

6.0 

14.697 

14.700 

8.518 

8.516 

3.449 

3.448 

4.349 

4.350 

-2.137 

-2.138 

5.757 

5.734 

Sources  of  the  values  of  thermodynamic  parameters: 

Thermal  Expansion  Coefficients:  R.  K.  Kirby,  J.  Res.  NBS,  A,  Phys.  Chem.  71A,  363  (1967). 
Specific  Heat  Value:  J.  S.  Arthur,  J.  Appl .  Phys.  21_,  732  (1950). 

Temperature  Derivatives  of  the  Compressibilities  at  1  atm,  Ref.  4.1. 


TABLE  4; 2 


Adiabatic  Elastic  Constants  (in  kilobars)  of  Single  Crystal  Rutile  with  Pressure  at  298°K  Obtained  by 

4  8 

Manghanani  '  (M)  and  by  Iterative  Procedure  (IP)  from  Manghanani's  Data. 


Pressure 

C11 

C33 

C44 

'  C66 

C12 

C13 

Kbar 

M 

Ip 

M 

IP 

M 

Ip 

M 

IP 

M 

IP 

M 

IP 

0.001 

2714.3 

2714.3 

4839.5 

4839.5 

1244.3 

1244.3 

1947.7 

1947.7 

1779.6 

1779.6 

1495.7 

1495.7 

1.25 

2722.4 

2722.4 

4849.9 

4850.8 

1245.7 

1245.2 

1955.7 

1955.7 

1791.0 

1791.0 

1502.0 

1503.4 

2.50 

2730.5 

2730.4 

4860.3 

4862.1 

1247.1 

1246.3 

1963.8 

1963.7 

1802.3 

1802.3 

1508.3 

1509.6 

5.00 

2746.7 

2746.6 

4881.2 

4884.7 

1249.8 

1248.4 

1979.9 

1979.7 

1825.1 

1825.1 

1520.8 

1522.2 

7.50 

2762.8 

2862.8 

4902.0 

4907.3 

1252.6 

1250.5 

1996.0 

1995.9 

1847.8 

1847.8 

1533.4 

1534.8 

Sources  of  the  values  of  thermodynamic  parameters: 

Thermal  Expansion  Coefficients:  R.  K.  Kirby,  M.  Res.  NBS,  A,  Phys.  Chem.  71A,  363  (1967). 
..  .  Specific  Heat  Value:  J.  S.  Arthur,  J.  Appl .  Phys.  21_,  732  (1950). 

Temperature  Derivatives  of  the  Compressibilities  at  1  atm.  Ref.  4.1. 
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elastic  constants  by  about  1%.  Section  IV  extends  the  above  iterative  procedure  toin- 
cTndp  data  and  constants  for  crystals  of  arbitrary  symmetry.  It  is  applied  to  calcite 


HD  form  i  A  70  and  rutite  for  illustration. 
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